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Abstract 

Using data collected by the CLEO III detector at CESR, we report on measurements of T(IS') 
decays to charmonium final states. The data sample used for this analysis consists of 21.2 x 10^ 
T(IS') decays, representing about 35 times more data than previous CLEO T(IS') data samples. 
We present substantially improved measurements of the branching fraction B{T{1S) J/ifj + X) 
using J/il) — > jJi^ and J/^jj ^ e^e^ decays. The branching fractions for these two modes 
are averaged, thereby obtaining: ^(T(15') J/^j + X) = (6.4 ± 0.4(stat)±0.6(syst)) x lO"^. A 
greatly improved measurement of the J/'^ momentum distribution is presented and indicates a 
spectrum which is much softer than predicted by the color-octet model and somewhat softer 
than the color-singlet model. First measurements of the J /ip polarization and production an- 
gle are also presented. In addition, we report on the first observation of T(IS') — > i^{2S) + X 
and evidence for T(IS') — > XcJ + X- Their branching fractions axe measured relative to 
B{T{IS) J/^ + X) and are found to be: ^fflSw/^S^ = ^-^l ± 0.11(stat)±0.08(syst), 
= 0.35±0.08(stat)±0.06(syst), Iffgff^J/^^^^) = 0.52±0.12(stat)±0.09(syst), and 

B(Y(^i5HJ/V'+$) ^ '^■^ confidence level. The resulting feed-down contributions to J/ if) are 

(24±6(stat)±5(syst))% for ^(25"), (ll±3(stat)±2(syst))% for Xci, (10±2(stat)±2(syst))% for Xc2, 
and < 8.2% at 90% confidence level for XcO- These measurements (apart from Xco) are about a 
factor of two larger than expected based on the color-octet model. 



*0n leave of absence from University of Chicago. 
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I. INTRODUCTION 



Charmonium has played a crucial role in the recent history of particle physics. It has been 
nearly 30 years since its discovery in both e^e~ interactions [1] and in collisions of protons 
on a beryllium target [2]. Over the last two decades, the charmonium and bottomonium 
systems have served as a laboratory for testing QCD. In the weak sector, charmonium also 
serves as a critical tool in extracting CKM [3] phases in 5-meson decays. However, even after 
30 years of studying cc systems, we still lack a complete understanding of their production 
mechanisms in glue-rich environments. 

About a decade ago, the CDF experiment reported production rates of charmonium in 
proton-antiproton collisions which exceeded the existing theoretical calculations by a factor 
of about 10 for J/iIj and about a factor of 50 for iIj{2S) [4]. An explanation of this excess was 
given by the so-called color-octet mechanism [5], whereby a gluon fragments into a color- 
octet cc pair, which then evolves non-perturbatively into a color-singlet by emission of 
a soft gluon. The size of this non-perturbative matrix element is not predicted and was 
determined by a fit to the CDF data. Because of the glue-rich environment and the a1 
suppression of the color-singlet process relative to the color-octet process, it was argued [5] 
that the latter contribution is likely to be important. While this model can explain the rate 
and momentum spectrum of J /if) and ip{2S) production at the Tcvatron it appears that it 
does not properly describe recent J/ip polarization data from CDF [6]. Furthermore, when 
the same matrix elements determined at CDF are applied to photoproduction of J/ip at 
HERA, the color-octet contribution is about a factor of ten too large [7]. 

Over the last several years, the role of the color-octet mechanism in J/ip production 
in e+e~ collisions has been under theoretical study [8]. The dynamics of the color-octet 
processes are expected to give rise to significant differences in the J/ ip momentum spectrum 
and production angle as compared to color-singlet production. Recently, both BaBaR [9] 
and Belle [10] have reported measurements of the cross section and the momentum spectra 
of J/ip^s in e+e" collisions on the T(45') {B decays excluded) or just below the T(45') 
{i.e., in the continuum). BaBaR measures a{e^e~ J/ip + X) = (2.52 ± 0.21 ± 0.21) 
pb, whereas Belle measures a number which is 40% lower, (1.47±0.10±0.13) pb (about 
3 standard deviations below than the BaBaR result). The two experiments both observe 
similar shapes for the J/ip momentum spectrum, which are softer than the predictions of 
the color-octet model [8] which predict a peaking of the J/ip momentum spectrum near 
the kinematic endpoint. However, recent theoretical studies of the color-octet subprocesses, 
e'^e~ J/ip + g [11] and e'^e" — > J/ip + gg [12], show that the perturbative expansion 
breaks down near the kinematic endpoint, and the authors appeal to soft-coUinear effective 
theory (SCET) to systematically include the non-perturbative effects. In Ref. [11], it is 
shown that by using SCET the color-octet model predictions can be sufficiently softened 
and reasonably good agreement with e+e~ — >• J/ip + X data can be achieved, although the 
calculation is not completely predictive because it uses a shape function which is fit to the 
e'^e~ — > J/ip + X data [9, 10]. Belle also reports on production of ip{2S) in e+e" collisions, 
with a measured ratio o-(e+e- ip{2S) + X)/o-(e+e- J/ipdhect + X) = 0.93 ± O.UtoH- 
[10]. That is, the production rates for J/ip and ^(25') in e"'"e~ collisions are approximately 
equal. The color-singlet mechanism can yield the ip{2S) final state, but the expected ratio is 
<10% [17]. Belle has extended their inclusive J/tp analysis to search for associated charmed 
particles, and they find: a{e+e- j/^jc!c)/a{e+e- J/^ + X) = 0.59t[J:l| ± 0.12 [13]. 
The color-octet contribution is expected to be at the level of 1% [14] of the inclusive rate. 
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The disagreement here indicates that the production mechanisms of charmonium are not 
well understood, and more theoretical and experimental input is required. 

Several theoretical papers [14, 15] have suggested that the study of J/%Ij production in 
T(15') decays could provide an alternate probe of the charmonium system in that the T(15') 
decay provides a glue-rich environment in which J /^ijj mesons can be produced abundantly 
through the color-octet mechanism. The kinematics of such J/ip^s are expected to exhibit 
signatures distinct from other production mechanisms, such as a peak in the J/ip momentum 
spectrum near the kinematic endpoint. The predicted branching fraction from color-octet 
processes is B{T(1S) — > J/t/j + X) — 6.2 x 10~^ [15], with approximately 10% feed-down 
expected from ip{2S) and another 10% from Xcj [16] (summed over all J). Color-singlet 
processes, such as T(IS') J/ip -\- gg start at a^, and are therefore suppressed relative to 
color-octet processes which enter at a^. However, computations of the color-singlet process 
T(IS') — > J/i/jccg+X [17] indicate a sizeable branching fraction of 5.9x10"^, with about 10% 
coming from ip{2S) feed-down. The enhancement here arises because the non-perturbative 
color-singlet matrix element for cc J/ip may be 210-360 times larger than the corre- 
sponding color-octet matrix element, which is enough to compensate for the perturbative 
suppression. Moreover, unlike the color-octet processes, this process inherently results in a 
soft J/ip momentum spectrum because of the two additional charm quarks in the final state. 
As a result, the J/ip momentum cannot exceed about 3.3 GeV/c in this process. Therefore, 
while the color-octet and color-singlet processes give similar predictions for total rate, their 
momentum distributions are significantly different. Figure 1 shows Feynman diagrams for 
(a) two of the more important color-octet processes and (b) the T(15') — > J/t/jccg + X color- 
singlet diagram. It should be noted that color-singlet production would also be signalled by 
the presence of additional charmed particles (open charm) in association with the J/ip. To 
capitalize on the small yield of J/ip^s in T{1S) decay, many D decay channels, both inclusive 
and exclusive will need to be explored. We therefore relegate the search for open charm in 
association with J/ip in T(IS') decay to a future report. 

The process T(IS') J/ip + X has been previously observed by CLEO [18], where the 
branching fraction was measured to be (l.l±0.4±0.2)x 10^'^ based on ^20 observed events. 
CLEO also reported a soft momentum spectrum for the J/ip, albeit with limited statistical 
precision. The ARGUS Collaboration reported an upper limit of 0.68x10"^ [19] at 90% 
confidence level. 

The CLEO Collaboration has collected large data samples on the T(nS') resonances and 
currently has the world's largest samples of T{1S), T{2S) and T(3S') decays. Consequently, 
CLEO is in a unique position to help clarify the roles of color-singlet and color-octet models 
in J ftp production. 

In this paper, we present vastly improved measurements of the rate, momentum spectrum, 
and angular distributions in T(IS') J /tp + X decays. We also present first observations 
of the decays T(15') i'{2S) + X and evidence for T(15') — > Xci,2 + X. The paper is 
organized as follows. In Section II we discuss the data samples used, the J/'4> backgrounds, 
event selection and J /ip reconstruction. Section III details the measurement of the T(IS') — >■ 
J/ip + X branching fraction and momentum spectrum in T(IS') decays. This section also 
includes a measurement of the cross section a{e'^e~ — > J/ip + X) using data on and below 
the T(4S') resonance, which is used to estimate and subtract the continuum contribution 
at the T(IS'). The section concludes with an examination of some event-level distributions. 
Section IV presents the measurement of the T(IS') ip{2S) + X branching fraction. The 
report then discusses in Section V the measurement B{T{1S) — > XcJ + For each of 



4 



1631004-098 



bb 





cc 



bb 




h-nnnnmp- 



( a ) Color - Octet Diagrams 



bb 




cc 





c^ 

c^ 
c 



cc 



( b ) Color - Singlet Diagram 

FIG. 1: Feynman diagrams for production of charmonium in T{1S) decays from (a) color-octet 
processes and (b) color-singlet T{1S) J/ipccg + X. For the color-octet processes, the J/V' is 
produced in a color-octet and becomes a color-singlet through emission of a soft gluon. 

these analyses, we present a cross-check by measuring the corresponding branching fraction 
in i?-meson decay. Lastly, we discuss in Section VI the systematic uncertainties in each of 
these analyses, and the paper is concluded in Section VII. 



II. DATA SAMPLES, BACKGROUNDS, EVENT SELECTION, AND Jf-ip RE- 
CONSTRUCTION 

The analysis presented here uses data collected using the CLEO III detector [20]. The 
primary data sample includes 1.2 fb^^ of data collected on the T(15'), and amounts to 
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(21.2±0.2)xl0^ T(IS') decays. For background determinations and systematic checks, we 
also utilize 5.0 fb^^ of data on the T(4S') resonance (10.4 million 5-meson decays) and 2.3 
fb~^ just below (~10.56 GeV) the T(45') resonance. We also use the on-T(45') data for 
cross-checks on charmonium yields in S-meson decays. 

The backgrounds to J/ip on the T(IS') are: (a) radiative Bhabha events, (b) 

77 fusion producing Xcj which subsequently produces J/ipj, (c) radiative return processes 
such as e+e" J/ipJ or e"'"e~ — >• ip{2S)'y, and (d) continuum production (e+e" —>■ J/ip + 
X). Various event selection requirements are targeted at reducing or eliminating these 
backgrounds. Radiative Bhabha events produce background in the J/ip mass region when 
one of the hard leptons is combined with a soft lepton from the converted photon. Such 
events are suppressed by requiring that the invariant mass of either electron from the J/ip ^ 
e^e~ candidate with any other electron in the event have Mee > 100 MeV/c^. Events 
produced through 77 — > XcJ fusion typically only leave two charged tracks in the CLEO 
III detector, and these events are therefore easily rejected by a requirement of at least 3 
charged tracks. The radiative return backgrounds are suppressed through event selection 
criteria which take advantage of the special kinematics of these processes, namely a low 
particle multiplicity coupled with either the detection of a high energy photon (Ri4 GeV) or 
large missing event momentum. Events are required to have their missing event momentum 
magnitude, Pcv < 3.75 GeV/c, or, if the number of charged tracks, A^trk < 4, we require 
Pcv < 2.0 GeV/c. When the high energy photon is detected (or an e+e^ pair with invariant 
mass less than 100 MeV/c^), the event is vetoed if iV^rk < 4 and the (converted) photon 
has energy greater than 3.75 GeV. The remaining background from these three sources 
to the T(IS') — s> J + X signal is negligible. However, because of the small signal in 
T(IS') 'ip{'2S) + X, the remaining background cannot be neglected. This background is 
determined using the evtgen MC followed by a GEANT-based detector simulation, and 
the resulting contribution is subtracted from the observed yields. 

Continuum background is reduced by requiring that the second Fox- Wolfram moment [21] , 
i?2 < 0.6. The remainder of this background is estimated using T(4S') data, and is statis- 
tically subtracted from the observed T(15') yields. The estimate of this background is 
discussed in Sections III B and III E. 

Candidate J/ip^s are formed by pairing oppositely charged electron or muon candidates. 
These charged track candidates are required to have momentum in the range from 0.1 to 
5.3 GeV/c and have at least 50% of the maximum number of expected hits in the tracking 
system. We also require these tracks to be consistent with coming from the interaction 
point in three dimensions. Electron candidates are additionally required to have a shower 
profile which is consistent with expectations for an electron and an energy deposition in 
the calorimeter, Ee, which is compatible with its measured momentum, pe, by requiring 
0.85 < Eg/pe < 1.15. For these electron candidates, we correct for radiated photons by 
adding back the momentum of the highest energy photon which lies within a 5° cone of the 
initial particle direction. Muon candidates are formed using charged tracks which penetrate 
at least 3 hadronic interaction lengths of iron absorber in the muon chambers [20] . 
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FIG. 2: Dilepton invariant mass distributions for (a) J /ip fJ-^f^^ and (b) J/ip e^e~ candidates 
for data taken on the T(IS') resonance (points) and data taken just below the T(4S') resonance 
(shaded). The T(4S) distributions are scaled to account for the different integrated luminosities 
and center of mass energies for the two data samples. 



III. MEASUREMENTS OF T(IS') ^ J/tp + X 
A. J/ip Mass Distributions in the T(IS') Data 

Figure 2 shows the invariant mass distribution of J/ip candidates for (a) J/tp ^ A*'*'/^" 
and (b) J/i(j — *• e~^e~ in the T(IS') on-resonance sample satisfying all selection criteria. 
The shaded histograms show the corresponding distributions for T{4S) continuum data, 
scaled by a factor of 0.65, which accounts for the differences in luminosities and center of 
mass energies. The mass distributions from the T(15') data set are fit to the sum of a 
linear background and a Gaussian signal whose means and widths are allowed to float. The 
fitted peaks have a resolution of 13.4 MeV/c^ and 14.2 MeV/c^ for the J/ip —>■ fi^fi~ and 
J/ip e^e~ channels, respectively. The fitted yields are 399±25 J /ip ^ /^^A*^ s-nd 449±27 
J/ip e^e~ signals events. 

To study the momentum distribution, we divide the data into bins of scaled momentum, 
x, where x = pj/^/pmax- Here, Pmax = (1/2a/s)(s — Mj^^) is the maximum J/ijj momentum 
assuming the J/ip is recoiling against a massless particle, s is the square of the center of 
mass energy, pj/^ is the momentum of the J/ip candidate and Mj/^ is the J/ip mass [22]. 
The data are binned in intervals of Ax = 0.2. This scaled momentum variable removes 
the beam-energy dependence which is useful in comparing spectra on the T(IS') and the 
T(4S'). The invariant mass distributions for J/ip — > and J/tp — > e+e~ for T(IS') data 

in bins of x are shown in Fig. 3 and Fig. 4, respectively. If the x distribution has a sharp 
peak near the kinematic endpoint, there may be smearing into the x > 1.0 region. The 
absence of any signal in the 1.0 < x < 1.2 bin shows that all events are containe within 
the physically-allowed region. A simulation of the J/ip signal (see Section IIIC) indicates 
that the widths of the invariant mass distributions are independent of J/ ip momentum, and 
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FIG. 3: Invariant mass distributions for J/'^ ^ jjL^ix candidates in x bins of size 0.2 for data 
taken on the T(IS') resonance. 



therefore these distributions are fit using a width fixed to the values obtained from the full 
sample. 



B. Candidate J/ if) Mass Distributions in the t{AS) Data 

The continuum contribution to the T(15') — > J/ip + X signal is estimated using data 
taken on and below the T(45'). This measurement is interesting in itself in light of the 

disagreement in the rates for e+e" J/ip + X measured by BaBaR [9] and Belle [10]. We 
employ the same event selection criteria as for the data taken on the T(IS'), except that 
for the on-T(45') data, we require the J/ip to have momentum larger than 2 GeV/c, which 
eliminates contributions from S-meson decay. 

The measured signal for e^e~ J /ip + X below the T(4S') is shown in Fig. 5 for (a) 
J/ip — > and (b) J /ip ^ e+e^. The fitted numbers of events arc 112±17 {J /ip — > Z^"*"/^") 

and 116±19 {J/ip e'^e"). The corresponding distributions for data taken on the T(45') 
resonance are shown in Fig. 6. The fitted yields are 130±17 J/tp ^ and 193±24 

J/'0 — e^e~ events. The yields per unit luminosity are statistically compatible, after 
correcting the on-T(45') yield for the 2 GeV/c momentum requirement. The correction is 
determined from the J/ip momentum spectrum from the below-T(45') continuum and is 
estimated to be (25 ± 6)%. 



C. J/^ Reconstruction Efficiency 

The data are corrected for geometric acceptance and analysis requirements using the 
PYTHIA Monte Carlo [23] and a GEANT-based detector simulation [24]. 
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FIG. 4: Invariant mass distributions for J/'^ ^ e^e" candidates in x bins of size 0.2 for data taken 
on the T(IS') resonance. 

The reconstruction efficiency as a function of x and cos ^j/^, where Ojj^ is the polar angle 
of the J /ip in the lab frame, is shown in Fig. 7. The circular points are for J ^ and 
the triangles are for J /ip ^ e^e" . The efficiencies decrease slightly with increasing momen- 
tum and I cos 9 J /^\ and have average values of (40±2)% for J/ip and (50±2)% for 
J/%1) — > e'^e~ . The small drop in efficiency with momentum is a result of not reconstructing 
the softer lepton which is emitted backward in the J / ip rest frame. The lower J/ip ^ fi~^fi~ 
reconstruction efficiency is due to the requirement that both muons penetrate at least three 
layers of iron absorber, which limits the muon momentum to be larger than about 1 GeV/c. 

The momentum distributions of both the T(IS') — > J/ip -\- X signal as well as the on- 
T(4S') and below-T(4yS') yields are corrected using these x-dependent efficiencies. This is 
justified since the reconstruction efficiency is not sensitive to small differences in the event 
environment (between T{\S) — >■ J/ip + X and e+e^ J/ip + X). Continuum-produced 
J/'^'s have a similar charged-track multiplicity to T(15') J/ip + X (^7-8) and the T(15') 
data peak at low R2 (see Section IIIH, and figures therein) as do R2 measurements in the 
continuum (see Fig.l (c),(d) in Ref. [9]). 

D. Corrected J/ip Momentum Distributions on and just below the T{AS) Reso- 
nance 

The resulting differential cross sections, da/dx, versus x, are shown in Fig. 8 using the 
combined on-T(45') and below-T(45') data. The circles represent J/ip ^ l^'^l-i'' a-^d the 
triangles are J/ip ^ e+e~. The distributions clearly peak at large x values with a mean 
of about 0.7. Integrating these distributions, and using B{J /ijj — > /+/^)=5.9% [22], we find 
cT(e+e- J/i) + X) = 2.0 ± 0.2(stat) pb for J/%} fi+fi' and a{e+e- J/tp + X) = 
1.7±0.2(stat) pb for J/ip ^ e+e~. Combining these results we obtain a{e'^e~' — > J/ip+X) — 
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FIG. 5: Invariant mass distribution for (a) J/ijj ^ /u^/x and (b) J/V' — e+e in data taken below 
the T{4S). The data are integrated over all momenta. 



1.9 ± 0.2(stat) pb. The results using the different lepton species are consistent with one 
another and lie between the BaBaR and Belle measurements of (2.52 ± 0.21 ± 0.21) pb and 
(1.47±0.10±0.13) pb, respectively. The rates found in the continuum are about a factor of 
6-7 lower than on the T(15'). 



E. Extrapolation of the T(4S) results to the T(15) 

The extrapolation of the differential cross section for e~^e~ J/i/j + X on and below the 
T{4S) (see Fig. 8) to the T(IS') requires that we take into account the differences between 
these two samples, and includes two factors (other than the luminosity scaling): the ratio 
of partonic cross sections for e'^e" — > J/ip + X and a phase space correction for producing 
the J /ip + X final state. For the former, we assume 1/s scaling, since the process proceeds 
through a virtual photon, and therefore the parton-level cross section at 9.46 GeV is 1.25 
times larger than at 10.58 GeV. For the phase space extrapolation, we bound this factor at 
unity by assuming the phase space at 9.46 GeV is equal to that at 10.58 GeV. To obtain 
a lower bound, we assume that the J/i/^'s are always produced in association with a pair 
of D mesons, which has a significantly reduced phase space at 9.46 GeV as compared to 
10.58 GeV. Using PYTHIA, we determine that the probability, of producing J/ipDD at 9.46 
GeV is 55% of the corresponding value at 10.58 GeV. Using these values as extremes, and 
assuming that the "true" value has a fiat probability of lying somewhere in that interval, 
we estimate the phase space ratio is 0.78 ± 0.13. Combining the two factors, we determine 
the continuum extrapolation factor, /cont = 0.98 ± 0.16. 

For the J/ijj momentum spectrum in T(15') decays, we are primarily interested in the 
shape for the gluonic intermediate states. The qq intermediate state, which proceeds through 
the coupling of the T(15') to a virtual photon, is assumed to have the same shape in x as in 
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FIG. 6: Invariant mass distributions for (a) J/ip ^ and (b) J/il^ 

the T(4S') resonance. To reject J/tp^s from B decay, we require the momentum of the J/ip to be 
larger than 2.0 GeV/c. 
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FIG. 7: Efficiency for reconstructing J/ip's in T{1S) decays as a function of (a) scaled J/xp mo- 
mentum, and (b) cosOj/^. The circles are for J/ijj — and the triangles are for J/ip — e+e~. 
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FIG. 8: Distributions in x for e^e^ — > J/^ + X using the combined data taken on and just below 
the T(45). The circles show the results obtained using J /if: ^ and the triangles show the 

corresponding distribution obtained using J/^ ^ e^e~ . 



e~^e~ J/ip + X, and therefore is more closely related to the predictions for J/ip production 
in the continuum. Therefore for the purposes of the momentum spectrum, we subtract the 
expected T(IS') -^'y*^qq^J/ip + X contribution. This contribution is included for the 
branching fraction measurement. Any potential interference between the continuum and 
the T(IS') ^ 7* ^ qq contributions is neglected. We express the T(IS') — 7* ^ qq rate 
relative to the corresponding rate for e~^e~ J/ip + X using, 

Here, <Jr(is)-^x is shorthand for cr(e"*'e^ T(IS')) x B(T{1S) X). The measured value 
for 0"x(i5)-+/i+/i- is 0.555±0.022 nb [25]. In that same reference, the theoretical value for 
ae+e-^fj.+fj,- at 9.46 GeV is estimated to be 1.12 nb [25]. A more recent estimate based on 
the FPAIR Monte Carlo (MC) simulation [26] gives a larger cross section of about 1.38 nb. 
Taking the average of these two cross sections as our central value and half their difference 
as the uncertainty, we obtain (Je+e-^^+;_t-=1.25±0.13 nb. We therefore estimate that the 
T(l^) ^ 7* ^ gg ^ J/ip + X contribution is (44±5)% of the e+e" J/^ + X at ^^5=9.46 
GeV. Adding this contribution to the continuum extrapolation factor, /cont? we obtain an 
overall extrapolation factor for the x spectrum of fx = 1.41 ± 0.18. 
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FIG. 9: Differential cross sections in x for data taken on the T(IS') (points) and data taken 
both on and just below the T(4S) (triangles). The latter have been scaled to account for the 
T{1S) ^* ^ qq contribution, as discussed in the text. The upper figure shows the results 
obtained using J/jp ^ /^^/^~ and the lower figure shows the corresponding distributions obtained 
using J/ip — > e+e~. 



F. Corrected J/tp Momentum Distributions and Branching Fractions in X{1S) 
Data 

Figure 9 shows the differential cross sections in x for (a) J/ip —>■ fj^^fi^ and (b) J/tp ^ 
e^e~ using data taken on the T(IS') (solid circles) and averaged results from the data taken 
on and below the T(4S') (triangles). The latter have been scaled as discussed above to 
include both the continuum and T(15') 7* — > qq contributions. The differential cross 
section (versus x) for T(IS') J/ip + X is given by the difference of these two distributions, 
and reflects only the contributions from gluonic intermediate states. The results are shown 
in Fig. 10 using J/ip ^ A*'*'/^" (circles) and J/ip e^e^ (triangles). The figure also 
shows the theoretical predictions of the color-octet [14] (solid line) and the color-singlet 
T(IS') J /ip + ccg [17] (dashed line) model. 

The branching fraction for T(IS') J/tp + X is computed by integrating the differential 
cross section distribution. We only subtract the expected continuum contribution (we use 
/cont as our extrapolation factor as opposed to /^.), so that the branching fraction includes 
the three intermediate hadronic states: ggg, 991^*^ and qq. The resulting branching fractions 
for the J/ip ^ and J/ip ^ e^e~ final states are: 



i3^^(T(l^) J/i) + X) = (6.9 ± 0.5(stat)) x 10"^ 

BeeiyilS) J/i) + X) = (6.1 ± 0.5(stat)) x 10"^ (2) 

Systematic uncertainties are presented in Section VI. Using rtot(T'(15')) = 53.0 ± 
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FIG. 10: Differential cross sections in x for T(IS') J/ip+X obtained using J/ip ^ /^^/^^ (circles) 
and J/tl^ — > e'^e~ (triangles). We also show the theoretical expectations based on the color-octet 
(solid line) [14] and color-singlet (dashed line) [17] models. 

1.5) keV [22] , our measurement corresponds to partial widths, T ggg^ggry_i_qg of (36.6 ±2.8) eV 
and (32.3 ± 2.8) eV for the J/ip ^ A*'*'/^" and J/ip —>■ e^e~ channels, respectively. 

Subtracting the expected T(IS') ^ 7* — * gg contribution, we obtain ^ggg+gg^yM of (33.9 ± 
2.8) eV and (30.2 ± 2.3) eV. In other words, about 90% of the J/ip rate comes from the ggg 
and gg'y'^*^ intermediate states. The gg'j^*^ contribution is only expected to be at the level 
of about 5% [14] of the ggg rate. 

Theoretical estimates of this rate based only on color-octet contributions, which neglect 
the qq intermediate state, give a total branching fraction of 6.2 x 10~^ [14, 15]. Those 
predictions are in good agreement with the measurements reported here. On the other 
hand, our measured momentum spectrum is significantly softer than predicted by the color- 
octet model, which is expected to peak near the kinematic limit (see Fig. 10). However, 
it has been recently pointed out [11] that in a similar process, e^e^ J/ip + X, the 
non-relativistic calculations break down near the kinematic endpoint where there are large 
perturbative and non-perturbative corrections. These effects may be systematically treated 
using so-called Soft-CoUinear Effecive Theory (SCET) and are expected to soften the J/ip 
momentum spectrum. Using SCET, the shape of the measured J/ip momentum spectrum 
in e~^e~ — J/ip + X, which peaks near x ^ 0.7, was shown to be reproducible [11]. It will 
be interesting to see if these corrections, when applied to T(IS') — >• J/ip + X, can soften the 
color-octet predictions sufficiently to bring them into agreement with our data. 

Our measured rate is also consistent with the predictions of the color-singlet process 
T(IS') —>■ J / tpccg+X, which predicts a branching fraction of 5.9 x 10~^ and a soft momentum 
spectrum which peaks at x ~ 0.5 and has a kinematic limit of x < 0.9. While the data is 
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somewhat softer than the color-singlet predictions, it should be noted that this is a parton- 
level calculation and neglects the hadronization process. Inclusion of the hadronization 
of the charm quarks to charm hadrons softens the J/ip momentum spectrum, with more 
softenining occuring as the mass of the recoiling system increases. Further softening of the 
J/ip momentum spectrum occurs when including the feed-down of iIj{2S) and XcJ to J/ijj. 
Using a pythia simulation of the color-singlet process, we are able to obtain reasonably 
good agreement in the x < 0.6 region using our measured values for the feed-down from 
%1){2S), XcJ to J/ if] along with a reasonable, but arbitrary admixture of recoiling L>, and 
D** states. This is not necessarily evidence for color-singlet production, but it is suggestive. 

G. J/'0 Angular Distributions in T{1S) Data 

Angular distributions have the potential to differentiate the mechanisms for J/ip pro- 
duction in e^e~ collisions. Theoretical predictions for the production and helicity angle 
distributions for continuum production are available [8, 27, 28], but the calculations are yet 
to be done for T(15') decay. 

In the same spirit, we present distributions of the (polar) production angle, cosOj/^, of 
the J/ip and the helicity angle, cos^^hch where 6'hei is the angle between the positive lepton 
momentum in the J/ip rest frame and the J/ijj momentum in the lab frame. The efficiency- 
corrected J/ip ^ IJ'^l-i'' and J/ip ^ e'^e" channels are combined and shown in Fig. 11. Here, 
we subtract the expected T(IS') ^ qq and continuum contributions to extract the ggg 

and gg^^*^ shapes. The normalizations are arbitrary. The top figure shows the distribution of 
cos 6'hei for T(IS') data (points), pythia simulation (dashed), and a fit (dotted) to the form 
(1 + Acos^^hei), from which we find A = -0.48 ± 0.16 (xVdof=0.50). The bottom figure 
shows the distribution in cosOj/^ for T(IS') data (points), PYTHIA simulation (dashed) and a 
fit (dotted) to the form (1 + cos^ ^j/^), from which we find B = 0.01 ±0.16 (xVdof=1.44). 
The functional forms are the same as those used to describe the angular distributions for 
continuum production of J/ip mesons [8, 28]. The negative value of A indicates that the 
J/ip has a significant longitudinal polarization component (a positive value would indicate 
transverse polarization). For continuum production of J/i/j, the color-octet and color-singlet 
models differ greatly on their expectations for B at large values of scaled momentum, with 
B ~ —0.85 for the color-singlet model [28] and B ?a for the color-octet model [8]. If a 
large difference persists for T(IS') decay, the production angle distribution could be useful 
in differentiating these two mechanisms. We note that the PYTHIA simulation (using default 
parameters), which produces J/ip via T(15') J/ip + DD, appears to be in reasonable 
agreement with data. 

H. Event-Level Distributions 

Additional information on the T(15') — > J/t/j + X process can be obtained by studying 
various event-level distributions. We present distributions of: 

• Number of reconstructed charged tracks, A^trk (Fig- 12(a)) 

• Reconstructed neutral energy in the crystal calorimeter, -Eneu (Fig- 12(b)) 

• The second Fox-Wolfram moment, R2 (Fig. 12(c)) 
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FIG. 11: The helicity angular distributions, cos^hci (top) and production angle, cos^j/^ (bottom) 
of the J/V' in T(15) — > J /ijj+X. In each case, the points are the T(15) data, the dashed histogram 
is a PYTHIA simulation of T(15) J/^ + X and the dotted line is a fit to the T(IS') data as 
described in the text. 



• Invariant mass recoiling against the J/ip^ Mrecoil (Fig- 12(d)) 

In each case, we have performed a sideband subtraction, where the signal region is defined 
to be from 3.04 < Mi+i- < 3.14 GeV/c^ and the sideband region includes the M;+;- regions 
from 2.90 — 2.95 and 3.20 — 3.25 GeV/c^. The relatively small continuum contribution has 
not been subtracted. For each distribution, we also show the corresponding distribution 
from the PYTHIA MC simulation, which primarily produces a final state which consists 
of J/ipDD. The data are shown as points {J/ip — >■ fJ''^IJ'~ are circles and J/ip —>■ e^e~ are 
triangles) and the simulation is the histogram. We find that the charged particle multiplicity, 
which includes all charged particles, has a mean of about nine. The neutral energy, which 
comprises all energy in the calorimeter which is not associated with charged tracks, has an 
average of about (1.5 — 2.0) GeV, with most of the events having less than about 3.5 GeV. 
The Fox- Wolfram moment, i?2, peaks at low R2 in T(IS') — > J /ip + X data which indicates 
that these events tend to be more spherical than collimated (jetlike). The recoil mass can 
be used to discern whether there is another particle recoiling against the J/ip. It is defined 
by: 

Mrecoil = ^{^s-Eji^Y-p]j^, (3) 

where s is the square of the center-of-mass energy, and pj/^ and Ej/^ are the momentum 
and energy of the J / candidate. We do not observe any significant peaks in the recoil mass 
spectrum, indicating that the J/ip is usually not accompanied only by a second (bound) cc 
meson. The color-octet model predicts ~1% contribution to the inclusive rate whereas the 
color-singlet model does not predict the fraction of recoiling charm which is in the form of 
charmonium. 
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FIG. 12: Sideband-subtracted distributions of (a) number of charged tracks, A'trk; (b) neutral 
energy, -EneU) (c) the second Fox- Wolfram moment, i?2, and (d) the recoil mass, Mrecoil, for 
T(15) data and pythia simulation. The points (triangles) correspond to J/V' — > (•^/V' ~^ 

e^e~), and the histograms are the corresponding distributions from the simulation. 



I. Cross-check using B — > J/xp -\- X 

As a cross-check of our detector simulation and analysis procedure, we use the same 
tools to measure B{B J/ip + X) in T(45') data. The efficiencies for reconstructing J/ip 
in B —>■ J/ip + X events are about 5% lower than in T(IS') J/ip + X. In addition to 
the selection requirements described in Section II, we require the J/ip momentum to be 
less than 2.0 GeV/c. The yields are corrected for the expected continuum contribution, 
which is typically at the level of 1-2% of the B ^ J/iJj + X yield. The resulting branching 
fractions are found to be (1.17±0.03(stat))% and (1.14±0.02(stat))% for J/ip ^+ and 
J/ip ^ e^e~, respectively. These results are slightly higher than the world average value 
of (1.090 ± 0.035)% [22]. This difference is included as a systematic uncertainty in the J /ip 
reconstruction efficiency. 



IV. MEASUREMENTS OF T(IS') tp{2S) + X 
A. Measurements in T(IS') Data 

We search for T(15') ip{2S)+X using the decay mode ip{2S) J/ipTT^7i~ . Pion candi- 
dates must pass the previously mentioned track selection criteria and must have a measured 
energy loss in the tracking chambers within four standard deviations of the expected value. 
Using all pairs of oppositely charged pion candidates, we compute the invariant mass differ- 
ence, M(/+/^7r+7r^) — Af(/+/^), a quantity which has better resolution than M {J / ipn^ -k^) . 
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FIG. 13: Invariant mass difference M(l^l tt^tt ) — M{l^l ) for botli J/ip 
e+e~ candidates with invariant mass in the range from 3.00-3.14 GeV/c^. 
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We also require M{l^lr) to be in the range from 3.00 - 3.14 GeV/c^. The resulting dis- 
tribution for M(/^/^7r^7r~) — M{l^l^) is shown in Fig. 13, where we have summed over 
both Icpton species. The distribution is fit to the sum of a Gaussian signal shape and a 
second-order polynomial background. The width of the Gaussian is fixed to 2.5 MeV/c^, 
the value determined from B — > + X data. The 0.3 MeV intrinsic width [22] of the 

11^(23) is neghgible compared to the detector resolution, and is therefore ignored. The fitted 
yield is 56±11 events. The significance of the signal, S/ \/S + B, where S is the fitted signal 
and B is the estimated background, varies from 6 — 7, depending on whether B is estimated 
from the sidebands or the background function. 

The radiative return background, e+e~ — > '0(2S')7, is estimated using the evtgen [29] 
simulation package and published cross sections in Ref . [30] . The events are processed using 
GEANT and analyzed using the same analysis tools as the T(IS') data. The efficiency for these 
events to pass a loose hadronic event selection is (1.4±0.1)% for J/ip ^ l^^l^" decays and 
(8.4±0.6)% for J/il) e+e- decays. For these subsamples, a fraction, flf'"^^^^ = 0.40±0.09 
of J/V' — > A*"*"/^" decays also pass the analysis-specific selection criteria discussed in Section II. 
The corresponding fraction for J /if) e^e~ decays is y™dret,pass _ gj^g _|_ g g2. Xhe larger 
efficiency for the electron channel to pass the loose hadronic event selection results from the 
use of the calorimeter in defining this subsample of events. With the assumption that all 
data events which fail the analysis requirements are radiative return (discussed below) , the 
expected background contribution in data from radiative return is computed using 

((■radret,pass \ 
^TFiS^ Nfr^'K (4) 
hi / MC 

where the quantity in parentheses is the ratio of simulated radiative return events which 
pass the analysis-specific selection to those that are rejected. The quantity N^^^^'^^^ is the 
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FIG. 14: Event-level distributions for events rejected in the ip{2S) analysis. Distributions shown 
are (a) number of charged tracks, (b) neutral energy, (c) missing event momentum, and (d) cosine 
of the ■0(25') production angle. Solid points are rejected T(IS') — > il>{2S) +X candidate events and 
the histogram is the e^e~ i^i^S) + X radiative return simulation. 



number of rejected events in T(IS') data for each lepton species. We find jsf^^^^'^'^i = 5 
and N^^^^'^^i = 39 in the il^{2S) signal region, obtained through sideband subtraction. We 
therefore estimate radiative return contributions of 2.0±1.0 events and 5.9±1.0 events in 
the J/ip ^ and J/ip —>■ e~^e~ channels, respectively, and therefore a total of 7.9ibl.4 

background events from this source. 

The assumption that the rejected events in data are from radiative return is supported 
by comparing event-level distributions of these rejected events between data and simulated 
e~^e~ — ^ ?/'(2S')7 radiative return events. Figure 14 shows comparisons of (a) number of 
reconstructed charged tracks, A^trk, (b) neutral energy in the calorimeter, -Eneu? (c) missing 
event momentum, Pevent, and (d) the cosine of the angle between the i'{2S) direction and the 
beam axis, cos ^^(25). In all cases, the radiative return simulation reproduces the rejected 
events in T(IS') data, indicating that the rejected data events are mostly from radiative 
return. 

In Fig. 15 we show the analogous distributions for T(IS') data (points) passing all analysis 
selection requirements. The corresponding distributions from a MC simulation of %p{2S)DD 
are overlaid (histogram). These distributions are clearly quite different than the distributions 
for rejected events (see Fig. 14). 

The continuum background contribution is estimated using the measured cross section, 
(T(e+e- ^ ij{2S) +X) = 0.67 ± 0.09;[i:??, by Belle [10]. The expected number of ip{2S) 
continuum background events is then given by 
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FIG. 15: Event-level distributions for events accepted in the tp{2S) analysis. Distributions shown 
are (a) number of charged tracks, (b) neutral energy, (c) missing event momentum, and (d) cosine 
of the ip{2S) production angle. Solid points are accepted T(IS') — > ip{2S) + X candidate events 
and the histogram is a T(IS') — > 'ijj{2S)DD MC simulation. 



(5) 

We use the integrated luminosity C = 1.2 fb~\ and branching ratios of B{ip{2S) 
J/iIjti+ti-) = 0.318 ± 0.010 [22] and BiJ/ip l+l') = 0.059 ± 0.001. The reconstruction 
efficiencies are determined using a pythia simulation of T(IS') — >■ ip{2S) + X which is used 
to model the continuum as well as the signal, for the same reasons as mentioned previously. 
The efficiencies for both J/ip ^ ^~ and J/ip ^ e^e~ final states are nearly independent 

of momentum with average values, ef ^^^ = (17±1)% for / = /i and (23 ± 1)% for / = e. For 
the background extrapolation, we have assumed the same phase space suppression for ip{2S) 
as J/ip, and assign a 50% uncertainty to its value. We therefore expect a continuum back- 
ground contribution of 2.5±1.3 J/ip ^ f^^f^~ ci-nd 3.4±1.8 J/ifj ^ e+e~ events, which sum 
to 5.9±2.2 events. The error is dominated by the uncertainties in a{e^e~ — >■ ip{2S) + X) 
and /cont- As a consistency check, we have searched our 2.3 fb~^ continuum data sample for 
ip{2S), and we find 2.6t2!6 ^^cl 12±4 events in the /z/i and ee channels, respectively. Using 
the Belle cross section measurement, we would have expected 4.8 ± 0.7 /i/i and 6.5 ± 1.0 ee 
events, which is consistent with our observations. 

Combining the radiative return and continuum backgrounds, we estimate a total back- 
ground of 13.8±2.6 events. The uncertainty in the central value is included as a systematic 
uncertainty (see Section VI). 

We now compute B{T{1S) ^(25) + X). To reduce systematic uncertainty, the ip{2S) 
branching fraction is computed relative to i3(T(15') ^ J /ip + X) and is given by 
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BiTjlS) ^ ^(2g) + X) _ 1 / iVg - <S 




(6) 



where N^^^^f^ i^uitc) total number of il^(2S) {J/ip) signal candidates for lepton species 

; = e and Z = ^",t,E (^zfitck) is the expected 7/^(2^) ( J/^) background, and ejf^^^^ (e^f 
is the average reconstruction efficiency for ip{2S) (J/ip). A summary of the inputs used 
for the T(IS') i^i'^S) + X branching fraction computation are presented in column 4 of 
Table I. The table also shows in columns 2 and 3 the values for the J/ip ^ and 
J/%1) — > e^e~ channels separately. The event yields are consistent with one another. 



TABLE I: Various quantities relevant to the tp{2S) analysis as discussed in the text. Quantities 

included are the number of reconstructed candidates J/tp and ip{2S) candidates, their expected 
backgrounds, and efficiencies. The bottom two lines give the computed ratio of branching fractions, 
as discussed in the text. . 



Quantity 




ee 


Combined (ee + fj,fx) 




399±25 


449±27 


848±37 




53±11 


66±13 


119±17 


J/^ 
HI 


(40±2)% 


(50±2)% 


(45±2)% 




21±7 


35±8 


56±11 




4.5±1.6 


9.3±2.0 


13.8±2.6 




(17±1)% 


(23±1)% 


(20±1)% 


B{r{iS)^i>{2S)+x) 

e(T(15)^J/V+X) 


0.35±0.15 


0.46±0.15 


0.41±0.11 


B{r{lS)^i,\2S)+X)B(i>i2S)^J/i^+X) 


0.20±0.09 


0.27±0.09 


0.24±0.06 


B{r(is)^j/i>+x) 



In Table I, the number of J/ip background events in the T(IS') data is computed using 
the average measured cross section for e+e" J/ip + X of (1.9 ± 0.2) pb, the average 
efficiencies (also shown in Table I) and the continuum extrapolation factor /cont discussed 
in Section III E. The ratio of branching fractions is computed to be: 

BiTjlS) ^^{2S)+X) _ 

B{T{lS)^J/^ + X) -U-4i±U.iiistatj. (0 

That is, we find that the rate for T{IS) ^p{2S) +X is (41±11)% of the rate for T{IS) 
J/ip + X (systematic uncertainties are discussed in Section VI). It is interesting to note 
that in the T(45') continuum. Belle finds this ratio to be about 0.45 [10] with about a 20% 
relative uncertainty. Using B{ip{2S) J ji) + X) = (57.9±1.9)% [22], we find the feed-down 
contribution of ip{2S) to T(15') —>■ J/ip + X to be: 

B(T(1^) ^ i;{2S) + X)BmS) J/^ + X) ... 

This ratio is significantly higher than the expectations of either the color-octet model [14] 
or the color-singlet model in Ref. [17], each which predict a feed-down rate to be about 10%. 
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B. Cross-check by measuring B tJ^i^S) + X in T{4S) Data 

As a cross-check on our analysis, wc measure the yield for B — > ijj{2S) + X using 10.4 
million S-meson decays from the T(4iS') data sample, and use the same simulation tools to 
translate this into a branching fraction. The analysis techniques are also identical, except 
that we additionally require the momentum of the ip{2S) to be less than 1.5 GeV/c, which 
is the kinematic limit for its production in 5-meson decay. We find 129±16 and 144±18 
signal events in the J/ip ^ and J/ip —>■ e^e~ channels, respectively, of which 3±1 and 

4±1 events are expected from continuum background. The efficiencies are determined using 
generated PYTHIA hb events followed by a full detector simulation, and they are found to 
be (18±2)% for the J/ip ^ fi^fi~ channel and (24±2)% for the J/ip ^ e~^e~ channel. The 
branching fractions are measured to be (3.6±0.4(stat)) x 10~^ for the J/ijj — >■ iJi^ n~ channel 
and (3.0 ± 0.4(stat)) x 10~^ for the J/'4> ^ e+e~ channel. Thus we obtain good agreement 
with the world average value of (3.10 ± 0.24) x 10"^ [22]. 

V. MEASUREMENTS OF XcJ IN T(IS') AND r{AS) DATA 
A. Measurement of B(T(15) ^ XcJ + ^) 

We search for T(IS') — > XcJ + -'^ by reconstructing the XcJ J/'4>1 decay. Photon 
candidates are required to have energy >100 MeV, not be matched to a charged track and 
have a shower shape consistent with that of a photon. We also require that the invariant mass 
of this photon with any other photon in the event is greater than 2.5 standard deviations 
away from the 7r° mass of 135 MeV/c^ [22]. Photon candidates passing these selection 
criteria are combined with J/ if) candidates to form a Xc candidate. As done previously, 
we compute mass differences, M{l~^l~'y) — M{l~^l~), for J/ip candidates which have a mass 
in the range 3.00 < M{l^l~) < 3.14 GeV/c^. From this distribution, wc subtract the 
analogous distribution obtained from the J/%Ij sidebands, here defined as candidates with 
2.88 < M{l+l-) < 2.95 GeV/c^ or 3.20 < M{l+l-) < 3.27 GeV/c^. As was done for the 
'0(2(5'), we combine and average the J/ijj — >■ ijT and J/ip ^ e~^e~ channels. The invariant 
mass difference distribution, Mi+i-^, — Mi+i-, is shown in Fig. 16. The solid histogram is 
the T(15') data and the shaded histogram is the T(45') continuum, scaled by the ratio 
of luminosities. The T(15') data is fit using three Gaussians on top of an exponential 
background. The Gaussian means are restricted to lie within 5 MeV/c^ of the world average 
values of the differences between the Xcj and J/ip masses (see Table II) [22] and the widths 
are constrained to the values found from simulation: 12.5 McV/c^, 10.5 McV/c^ and 10.3 
MeV/c^ for J = 0, 1, 2. The larger width for J = is a result of the 16.2 MeV intrinsic 
width which is included in the simulation. The fitted yields are 0±13, 52±12, and 47±11 
events for J = 0, 1, 2, respectively. The significance of each signal is given by S/^/B, where 
S is the signal yield, and B is the estimated background within ±3 standard deviations 
of the fitted mean using the exponential background function. The significances are found 
to be 3.9 and 4.1 for the J—1 and J— 2 states, respectively. The averaged efficiencies for 
the J/ip n+n' and J/i/j e+e" channels are (27±1)%, (30±1)% and (28±1)% for J 
= 0, 1, 2 states, respectively. The efficiency for reconstructing the J = state includes a 
~2% loss of signal events due to events in the tails of the Breit-Wigner. Branching fractions 
for T(15') — s> XcJ + ^ are computed relative to T(15') — > J/ijj + X and are tabulated in 
Table II. The last column shows the measured fraction of J/'^'s which come from XcJ feed- 
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FIG. 16: Difference of invariant masses, M(l^l^j) — M(l^l^), for M{l~^l~) in the range from 3.00 
to 3.14 GeV/cP for I = n,e combined. The sohd histogram represents the T{1S) data and the 
shaded histogram is the below-T(45) data scaled by the ratio of luminosities. The arrows indicate 
the mass differences corresponding to the three Xc states. 

down, which is (11±3)% for the J = 1 state and (10±2)% for the J = 2 state. We only 
obtain upper limits on the J = state. Theoretical estimates of this ratio are at the level 
of 10% for the sum of all three Xcj states [14, 16]. The rates we report here are higher than 
those expectations. 

TABLE II: Measurements of the branching fractions for T(IS') — > Xcj + ^ with the relevant 
inputs. The table includes, from left to right, the Xcj states, the world average mass difference 
M^^j — Mj/^, the branching fractions B{xc,j — *■ J/i^l), the observed yields, the reconstruction 
efficiencies, the computed branching fractions relative to T(IS') J + X, and the computed 
feed-down to J/V'- For the XcO we show 90% confidence level upper limits. 



Mode 






-^events 


Eff.(%) 


B(T{15)-^J/»/'+X) 


Feed-down 




(MeV/c2) 


(%) 


(/x/x+ee) 


(/i/x+ee) 


(/u/x+ee) 


to Jji) 


XcO 


318 


1.11±0.15 


0±13 


(27±1)% 


< 5.9 


<0.065 


Xcl 


414 


31.6±2.7 


52±12 


(30±1)% 


0.35±0.08 


0.11±0.03 


Xc2 


459 


20.2±2.0 


47±11 


(28±1)% 


0.52±0.12 


0.10±0.02 



B. Cross-check using B — >^ XcJ + -X^ in T(4S) Data 

As a consistency check, we measure the branching fraction for B(B Xcj + X) using 
10.4 million B decays from the T(4S') data sample. We restrict the XcJ to have momentum 
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less than 1.6 GeV/c, which is the kinematic hmit from 5-meson decay. A clear signal is 
foimd only for the J = 1 state, for which the fitted yield is 347±35 events. The continuum 
background is negligible, and is therefore neglected. The reconstruction efficiency is deter- 
mined using a B ^ Xcj + X MC simulation, and is found to be (28±1)%. The branching 
fraction for B{B Xci + -^) is thus found to be (3.3 ±0.3(stat)) x 10~^, which is consistent 
with the PDG value of (3.6 ± 0.3) x lO'^ [22]. The 90% confidence level upper limit on 
B{B — > Xc2 + X) is 1.5x10"^, which does not conflict with the measured branching fraction 
of (1.3±0.4)xl0-^ [22, 31]. 

VI. ESTIMATES OF SYSTEMATIC UNCERTAINTY 
A. Uncertainties in T(1<S) J/tj} + X 

The branching fractions for B — > J/i/j + X using J/^ — > jji^ iJi~ and J /ip ^ e^e~ were 
shown in Section III. I to be higher than the world average values by 7% for J/ip ^ f^^f^~ 
and 5% for J/ip ^ e'^e~, which is taken as the systematic uncertainty in reconstructing 
these decays. We ascribe an additional uncertainty due to our limited knowledge of the 
final state in T(IS') J/ip + X and its modeling. This additional uncertainty is taken to 
be half the difference in the reconstruction efficiency obtained from our T(15') — > J/ip + X 
simulation and that obtained using the B ^ J/ip + X simulation. This results in additional 
contributions of 4% for J/i/j fi^fi^ and 6% for J/ip ^ e~^e~ . We also include an additional 
5% uncertainty in each due to limited MC statistics. We therefore estimate 9% systematic 
uncertainty in the J/ip ^ l^^l^^ reconstruction efficiency and 8% for J/ip ^ e+e^. 

The uncertainty in the signal yield is estimated by floating the Gaussian widths used in 
fltting each x bin. We flnd that the signal yield changes by 3% for both the J/t/j ^ l^'^l^' 
and J/i/j e+e~ analyses. The systematic error due to uncertainty in the shape of the 
background was estimated by comparing a linear background shape with an exponential. 
The difference is found to be 2% for J/-^ — > and 1% for J/ip ^ e+e^. 

Systematic uncertainty in the background subtraction comes from lack of precise knowl- 
edge of the continuum cross section (e^e~ ^ J/ip + X) and the error in the extrapolation 
from the T(4S') energy to the T(IS') energy. The latter includes uncertainties in the ratio of 
luminosities and the branching fraction i3(T(15') —>■ Our measurement of the rate 

for e+e~ — > J/'4> -\- X is uncertain at the level of 10% (statistical uncertainty only) and the 
extrapolation factor of 1.41±0.18 is uncertain at the level of 13%. We therefore estimate 
that the overall background rate is uncertain at the level of 16%. This uncertainty is propa- 
gated to an error in the branching fraction by shifting the central value for the background 
(see Fig. 8) up and then down by one standard deviation, and in each case, computing the 
change in the branching fraction from the nominal value. The corresponding shift in the 
branching fraction for T(IS') J/ip + X is found to be 6% for both J/ip fi~^/j~ and 
J/i/j ^ e+e", which is taken as the systematic uncertainty due to the background subtrac- 
tion. The uncertainty on the background subtraction does not have a significant effect on 
the general shape of the momentum distribution since its contribution is only about 10% of 
the total. 

We also include the uncertainty in the number of T(IS') decays, which is estimated to 
be 1% based on the number of T(IS') events and the uncertainty in the off-to-on T(15') 
luminosity ratio. We also include a 2% relative uncertainty in B{J/il} The total 

uncertainty is therefore found to be 12% for the J/V' — > A*"'"/^" channel and 11% for J/'^ — > 
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e^e channel. 

The systematics are itemized and shown in columns 2 and 3 in Table III. 

B. Uncertainties in T(1S) V(2S') + X Analysis 

For the "ip{'2,S) (as well as the Xcj) analysis, many of the systematic uncertainties cancel 
since these measurements are reported as a ratio with respect to the T(15') — > J/V' + 
X branching fraction. The uncertainty in the ip{2S) reconstruction efficiency comes from 
limited MC statistics (5%) and an imperfect understanding of T(15') ipi'^S) + X events. 
The uncertainty from the latter is taken to be the half the difference between the efficiency 
obtained using our default T(15') — > '^{'2.3) + X simulation and the B — > '4>{2S) + X 
simulation. The two simulations agree to within 1% in absolute value, which translates 
into an additional 5% relative systematic uncertainty in the reconstruction efficiency for 
each channel. We also include a systematic uncertainty of 1% per track for each of the two 
pions in the decay ^(25*) J/ipT^^T:' (2%). We therefore estimate that the non-cancelling 
systematic uncertainty in the tl^{2S) reconstruction efficiency is 7%. 

Uncertainty in the signal yield is estimated by shifting the Gaussian width up and down by 
20% about the central vahie (2.5 MeV/c^) and taking half of the average deviation, which 
results in a 7% systematic uncertainty. Uncertainty in signal yield due to the assumed 
background shape was estimated by fitting the background to the alternate functional form: 
^4(1 — B cxp{—Cx)) (the default is a second-order polynomial). The yield differs by 6%, 
which is taken as the associated uncertainty. 

Uncertainty due to the background subtraction is estimated by considering a ±50% 
change in the expected background contribution, which is about two standard deviations. 
The resulting systematic uncertainty in the branching ratio is 15%. Uncertainty in the 
branching fraction B{iIj{2S) ipTc^n^) contributes 5%. We therefore estimate a total sys- 
tematic uncertainty of 20% in the ip{2S) branching fraction ratio. 

C. Uncertainties in T(15) XcJ + X Analyses 

The uncertainty in the efficiencies for reconstructing XcJ J/i^l is taken to be half 
the difference in the efficiencies for reconstructing Xcj in B decays at the T(4S') versus in 
T(IS') decays (8%). We include an additional uncertainty of 8% to reflect the lower value 
we obtain for B{B —>■ Xci + ^) as compared to the world average. This also accounts for 
any possible systematic uncertainty in the photon reconstruction efficiency. We attribute 
a 3% uncertainty for each due to limited MC statistics. We therefore estimate an overall 
systematic uncertainty in the reconstruction efficiency of the Xcj + X final state of 12%. 

The uncertainty in the signal yield is obtained by allowing the Gaussian widths to float, 
from which we obtain differences of 6% and 1% for the J = 1 and J = 2 states. The 
uncertainty from the background determination is estimated by using different ranges over 
which to fit the backgroTind. We find maximum variations of 12%, of which we take half 
as the associated systematic uncertainty (6%). Since there is no evidence of any Xcj signal 
in the continuum, we do not consider this as a source of systematic uncertainty. Lastly, we 
include uncertainties in the XcJ J/ "4^1 branching fractions [22] of 20%, 10% and 11% for 
J = 0, 1 and 2 states, respectively. We therefore estimate total systematic uncertainties of 
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18% for Xcl and 17% for Xc2- The uncertainty on the hmits for Xco ai'c estimated to be 25%, 
and the upper hmits are increased by this amount to reflect this systematic uncertainty. 
Systematics uncertainties are hsted and summarized in Table III. 



TABLE III: Sources of systematic uncertainty in the T(IS') J/ip + X, T{IS) i){2S) + X and 
T(IS') — XcJ + X analyses. 



Source 


Value (%) 




J/^p 


J/il) ^ e^e 


^(25) 


XcO 


Xcl 


Xc2 


Reconstruction efficiency 


9% 


8% 


7% 


12% 


12% 


12% 


Signal Yield 


3% 


3% 


7% 




6 


1 


Background Shape 


2% 


1% 


6% 




6 


6 


Background subtraction 


6% 


6% 


15% 








#T(15) decays 


1% 


1% 










Error in B{ J/ip l^l^) 


2% 


2% 










Error in B{iIj{2S) ipir^TT^) 






5% 








Error in B{xcO,i,2 ^ J/'^Pl) 








20% 


10% 


11% 


Total 


12 


11 


20 


25 


18 


17 



VII. SUMMARY 

We present vastly improved measurements of the rates for production of charmonium in 
T(IS') decays over previous measurements. We have measured both the branching fraction 

for T(IS') — > .//'(/' + X and the scaled momentum distribution, as well as distributions in the 
polar angle and helicity in T(15') decay. We also report on first observations of the decays 
T{1S) i){2S) + X and evidence for T{IS) Xci,2 + X. 

The branching fractions for T(15') — > Jj^) ^ X are measured in both the Jj^) ^ A* """A*" 
and J/ip ^ e+e^ channels. Their braching fractions, B{T{1S) J/i/j + X) are measured to 
be (6.9 ± 0.5(stat) ± 0.8(syst)) x 10"^ and (6.1 ± 0.5(stat) ± 0.7(syst)) x 10-^ respectively. 
The two are averaged to obtain: 

B{T{1S) J/ij + X) = (6.4 ± 0.4(stat) ± 0.6(syst)) x 10"^ (9) 

We also measure the branching fraction B{T{1S) ^(25*) + X) relative to i3(T(15') — > 
J/^Ij + X), and find 

m=i||)^.0.4..0.n,.a,.0.0S,s..). (10) 

This report also presents the first evidence of the decay T(15') — > XcJ + ^- The branching 
fractions for all measured modes are summarized in Table lY. 

The T(15') — > J/t/j+X branching fraction is consistent with predictions of both the color- 
octet mechanism for J/ip production in T(IS') decays [14, 15] and color-singlet production via 
T(15') J/ipccg + X [17], each which predict a branching fraction at the level of 6xl0~^. 
The observed scaled momentum spectrum is relatively soft, peaking around x ~ 0.3, in 
contrast to J/ip^s produced in the continuum, which peak at about 0.7. The peaking at low 
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TABLE IV: Summary of measurements of the branching fractions for T(IS') to charmonium final 
states. The ipC^^S) and XcJ branching fractions are expressed relative to B{T{1S) — J/il^ + X). 
The last column show the feed-down contributions to T(liS') — > J/il) + X. For the XcO we show 
upper limits at 90% confidence level. Where uncertainties are shown, the first is statistical and the 
second is systematic. 



Final state, / 


B{TilS) ^ f + X) 


Feed-down to J/i/j 




(6.4 ±0.4 ±0.6) X 10-4 






B{r{lS) -f + X)/B{T{1S) - J/v + X) 




XcO 
Xcl 
Xc2 


0.41 ±0.11 ±0.08 

<7.4 
0.35 ± 0.08 ± 0.06 
0.52 ± 0.12 ± 0.09 


0.24 ±0.06 ±0.05 

<0.082 
0.11 ±0.03 ±0.02 
0.10 ±0.02 ±0.02 



momentum is in sharp contrast to the prediction of the color-octet model which predicts a 
peaking of ,t near 1. It is possible that incorporation of final state interactions could improve 
this agreement as was shown for e+e" — >■ J/ip + X [11]. The observed spectrum is closer to, 
although softer than, the expectation of the color-singlet process [17], T(15') — > J /ipccg + X, 
which peaks near x ~0.5. When this parton-level calculation is simulated using pythia we 
are able to achieve satisfactory agreement in the region x < 0.6 when hadronization of the 
recoiling charm quarks into charm hadrons is included, and our measured feed-downs of 
il'{2S) and XcJ to J/i^ are incorporated. 

The observation of T(15') — > ■0(25') + X is the first to a cc final state other than J/t/j. 
The feed-down to J/ip constitutes (24±6±5)% of the inclusive rate for T(IS') — >• J/ip + X, 
which is significantly larger than expected in either the color-octet [14] or color-singlet 
model [17], each which predict a feed-down to J/ip at the level of 10%. Our measured rates 
for T(15') — s> XcJ + X yield feed-down contributions of (11±3±2)% for the J—1 state and 
(10±2±2)% for the J— 2 state, which is also larger than the expected contribution of about 
10%, summed over J=0, 1, and 2 [14, 16, 17]. 

These measurements can shed additional light on the role of the color-octet and color- 
singlet mechanisms in producing charmonium, not only in T(15') decays but also in e+e" 
and pp collisions. In this regard, it would be of great interest to determine whether the same 
softening mechanism applied to the color-octet prediction for e~^e~ — > J/ip + X [11] can 
account for the J/ip momentum spectrum in T(15') J/i/j + X. Moreover, computation 
of the angular distributions for the color-octet and color-singlet mechanisms may provide 
additional discrimination between these two processes. Prom an experimental perspective, 
the additional information on the roles of color-singlet versus color-octet mechanisms may 
be obtained by measuring the ratio a{pp J/iIjcc + X)/a{pp — > J/ip + X) at the Tevatron. 
The unexpectedly large value for a{e'^e~ —>■ J/ipcc + X)/a{e'^e~ — >• J/ip + X) reported by 
Belle [13] may point to a large rate in pp collisions as well. 
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